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We study the critical temperature T c of FSF trilayers (F is a ferromagnet, S is a singlet superconductor), 
where the triplet superconducting component is generated at noncollinear magnetizations of the F layers. An 
exact numerical method is employed to calculate T c as a function of the trilayer parameters, in particular, mu- 
tual orientation of magnetizations. Analytically, we consider limiting cases. Our results determine conditions 
which are necessary for existence of recently investigated odd triplet superconductivity in SF multilayers. 

PACS: 74.45.+C, 74.78.Fk, 75.70.Cn, 74.62.Yb 



A striking feature of the proximity effect between 
singlet superconductors and nonhomogeneous ferromag- 
nets is the possibility of generating the triplet supercon- 
ducting component [1,2]. Recently, it was shown that 
the triplet component also arises in the case of several 
homogeneous but differently oriented ferromagnets [3]. 
Physically, the generating of the triplet component in 
SF systems [1-3] is similar to the case of magnetic su- 
perconductors [4]. 

In Ref. [3], the Josephson effect was studied having 
in mind that the superconductivity in the system is not 
suppressed by the ferromagnets. However, this issue re- 
quires separate study. 

Although the SF proximity effect is rather well stud- 
ied, the influence of the mutual orientation of F lay- 
ers magnetizations (exchange fields) on T c of layered 
SF structures has been mostly considered basing on 
the cases of parallel (P) and antiparallel (AP) align- 
ment [5-10]. At the same time, those are the only cases 
when the triplet component is absent. 

A FSF trilayer with homogeneous but noncollinear 
magnetizations of the F layers is the simplest example 
of a layered structure in which the triplet component 
is generated. The triplet component (correlations be- 
tween quasiparticles with parallel spins) arises as a re- 
sult of interplay between the Andreev reflections at the 
two SF interfaces. This mechanism is similar to the one 
described in Ref. [2], with the difference that instead 
of local magnetic inhomogeneity we deal with magnetic 
inhomogeneity of the structure as a whole. 

The critical temperature of the noncollinear FSF 
system was studied in Ref. [11]. However, in that work 
the triplet component was not taken into account. Thus 
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Fig.l. FSF trilayer. The system is the same as in 
Ref. [11]. The thickness of the S layer is 2d s , of each 
F layer — dj. The center of the S layer corresponds 
to x — 0. The thick arrows in the F layers denote the 
exchange fields h lying in the (y, z) plane. The angle 
between the in-plane exchange fields is 2a. 

calculation of T c in the noncollinear FSF trilayer is still 
an open question. 

In this letter we study the critical temperature of 
a FSF trilayer at arbitrary angle between the in-plane 
magnetizations (see Fig.l), which makes it necessary to 
take the triplet component into account. We reduce the 
problem to the form, which allows to apply general nu- 
merical methods developed in Refs. [12,13]. This form 
also leads to some general conclusions about T c and al- 
lows analytical progress in limiting cases. 

1. General description. We consider the dirty 
limit, which is described by the Usadel equations. Near 
T c , the Usadel equations are linearized and contain only 
the anomalous Green function F [1]: 
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Here D is the diffusion constant (D s and Df for the S 
and F layers), w„ — 7rT(2n + 1) are the Matsubara fre- 
quencies, and c?3 is the third Pauli matrix. The function 
F is a matrix in the spin space. The and fn compo- 
nents describe the triplet superconducting correlations. 
In the P and AP cases it is sufficient to consider only 
the scalar equation for the singlet component f^. 

Equation (1) is written in the general case when both 
pair potential and exchange field are present. In our sys- 
tem, in the F layers the pair potential is absent, A = 0, 
while 

H = h (<?2 sin a + CT3 cos a) (2) 

at the exchange field h = h(0, sina, cosa). h is the 
exchange energy, and a describes the direction of the 
in-plane magnetization. 

In the S layer, the exchange energy is zero, while the 
pair potential obeys the self-consistency equation 
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where T cs is the critical temperature of the S material. 
In the case of a single S layer, A can be chosen real. 

The boundary conditions at the outer surfaces of the 
trilayer are 

dF f /dx = 0, (4) 
while at the SF interfaces 

£ s (dF s /dx) =i£f(dF f /dx), 7 = p s €s/pf€f, (5) 
±Znb(dF f /dx)=F s -F f , T b = R b A/p f £f. (6) 



Here £ s (/) = y/ 'D s ^/2ttT cs and p s ^ are the coherence 
lengths and the normal state resistivities of the S and 
F metals, R b is the total resistance of the SF bound- 
ary, and A is its area. The ± sign in the l.h.s. of Eq. 
(6) refers to the left and right SF interface, respectively. 
The above boundary conditions were derived for SN in- 
terfaces [14] (N is a normal metal); their use in the SF 
case is justified by the small parameter h/Ep <C 1 (Ep 
is the Fermi energy). 

We expand the Green function F in the basis of the 
Pauli matrices <7j, i — 1,2,3, and the unity matrix do- 
It can be shown that the solution has the form 



F = f a + /1CT1 + /3CT3. 



The /o component is imaginary, while /1 and are real. 
The relations / (-w„) = -/o(w n ), /1 (-<*>„) = -/i(w„), 
f 3 (—uj n ) = fz(uj n ) make it sufficient to consider only 
positive Matsubara frequencies. 

The /1 component describes a special type of triplet 
condensate [1,3], odd in frequency [/i(— u n ) = — fi(co n )] 



and even in momentum, which is similar to the one pro- 
posed by Berezinskii [15]. It is independence on the 
momentum direction that allows the triplet condensate 
to survive in the diffusive limit. 

Equation (1) yields three coupled scalar equations 
(we consider u n > 0): 
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w n /o - ihf 3 cos a = 0, 



w„/i + hf 3 sma = 0, 



(8) 



Dd 2 h 
2 dx 2 



^ — ihfo cos a — hf\ sin a + A = 0. 



Analyzing symmetries implied by Eqs. (8) and geom- 
etry of the system, we conclude that fo(x) = fo(—x), 
fi(x) = -fi(-x), f 3 (x) = h{-x). Thus we can con- 
sider only one half of the system, say x < 0, while the 
boundary conditions at x — are 

df a /dx = 0, / 1 = 0, df 3 /dx = 0. (9) 

Below we shall use the following wave vectors: 

k f = J2uj n /D f , 



kh 



h/D f , 



k 



h = ^kj 



+ 2ik 2 h , k s = y / 2uj n /D s 



(10) 



The solution in the left F layer, satisfying the boundary 
condition (4), has the form 

Ff — C\ (zt?o sin a + u\ cos a) cosh [kf (x + d s + df)] + 



+C2 (<7o cos a. + iS\ sin a + ct 3 ) cosh kh (x + d. , + df) 



+ 



+C3 (cto cos a + id\ sin a — a 3) cosh k^ (x + d s + df) . 

(11) 

The matrix boundary condition (6) yields three scalar 
equations, which allow to express the coefficients C\, 
C2, C3 in terms of the components /o, /1, f 3 of the 
Green function on the S side of the FS interface: 

Ci = (-i/osina + /icosa)/(l+7(,-Af), 

C 2 = (/o cos a - ih sin a + / 3 ) /2 (1 + lb A h ) , (12) 

C 3 = (/o cos a - if 1 sin a - / 3 ) /2 (1 + -y b A* h ) , 



(7) where we have introduced the following notations: 



Af = kf£f tanh(fc/d/), A h = k h ^ f t&nh(k h d f ), 
Vf=jAf/(l + lb A f ), V h =jA h /(l + lb A h ). (13) 

Then the boundary condition (5) yields three scalar 
equations which entangle / , /1, and / 3 . Thus the Green 
function of the F layer is eliminated, and we obtain 
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equations for the S layer only. Moreover, we can proceed 
further, because in the S layer the unknown function 
A (a:) only enters the equation for the fa component [see 
Eqs. (8)]. At the same time, taking boundary condi- 
tions (9) into account, we can write /o = Bq cosh(fc s x), 
fi = Bi sinh(fc s x). Excluding B and B 1 , we arrive at 
the effective boundary condition for / 3 : 



Udh/dx) = Wf 3 , 



where 



W = ReV h + 



(lmV h y 



(14) 



(15) 



k s £ s A(a) +ReVh 
and the angular dependence is determined by 

fc s £ s ta,nh(k s d s ) + Vf [sin 2 a + tanh 2 (fc s <i s ) cos 2 a 



A = 



ksis [ c °s 2 a + tanh 2 (fc s (i s ) sin 2 a] + Vf ta,nh(k s d s ) 



Effectively, we obtain the following problem: 



w„>0 



(16) 



(17) 



D s d 2 .f 3 
2 dx 2 



UJnh + A = 0, 



dfo{-ds 
dx 



W(oj n )f 3 (-d s ), 



#3(0) 

dx 



(18) 
(19) 



- this is exactly the problem that was solved in Refs. 
[12, 13]. Inserting the new function W, we can use the 
methods developed in those works. At a = 0, Eq. (15) 
reproduces W from Refs. [12,13]. 

All information about the F layers is contained in 
a single function W, all information about the misori- 
entation angle — in its part A(a). Knowledge of W 
is already sufficient to draw several general conclusions 
about the behavior of T c . First, if the S layer is thick, 
i.e. d s S> £ s , then tanh(/c s d s ) w 1 at characteristic fre- 
quencies, and T c does not depend on a. Qualitatively, 
this happens because the effect of mutual orientation of 
the F layers is due to "interaction" between the two SF 
interfaces, which is efficient only in the case of thin S 
layer. Second, T c does not depend on df if df » £/. 
Qualitatively, this is due to the fact that the supercon- 
ducting correlations penetrate from the S to F layer only 
on the scale £/. 

The triplet component is "nonmonotonic" as a func- 
tion of a: it vanishes at a — and a = n/2 (P and 
AP case, respectively), and arises only between the two 
boundary values. However, the T c {a) dependence is al- 
ways monotonic. It can be directly proven from the 
monotonic behavior of A(a), and, hence, W . This rig- 
orously derived conclusion disproves the result obtained 
by the approximate single-mode method in Ref. [7], 
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Fig.2. Critical temperature T c vs. thickness of the F 
layers df, which is normalized on the wavelength of the 
singlet component oscillations = 2-njkh- Parame- 
ters d s li s = 1.2, h/nT cs = 6.8, 7 = 0.15, j b = 0.02 
correspond to Ref. [13]. The curves are calculated at 
different angles 2a between the in-plane exchange fields 
in the F layers. 

where it was claimed that T c in the AP configuration 
can be smaller than in the P case. 

Numerical results obtained by the methods devel- 
oped in Refs. [12,13], are shown in Figs. 2, 3. A question 
arises: why is there pronounced angular dependence in 
the case d s > £ s , when the S layer is not thin? The 
answer is that the condition d s <C £ s = \J D s /2irT cs is 
a sufficient condition of thin S layer, whereas the nec- 
essary condition is weaker: d s -C £ = y // D s /27rT c , since 
the characteristic energy for a particular system is 7rT c 
with its own value of T c . The two conditions become 
essentially different if T c is notably suppressed, and in 
this case T c can exhibit pronounced angular dependence 
at d s <§C £, while it is possible to have d s > £ s . 

Experimentally, the conditions for observing the an- 
gular dependence of T c are more easily met when T c 
is essentially (but not completely) suppressed. Accord- 
ingly, the effect of a on T c (df) dependence is most pro- 
nounced near the reentrant behavior. Experimental de- 
tection of such behavior was reported in Ref. [16]. 

2. Thin S layer. If d s <C £ s , then A is constant. 
The Usadel equation (18) can be solved, and the equa- 
tion determining T c takes the form 



T 
ln^ 



2:t7, VI — 



UJ„>0 



1 



U) r , 



WirT cs Z s /d s 



(20) 



where W is given by Eq. (15) with simplified function 
A(a): 

,2 



A = 



k 2 ^ s d s + Vf [sin 2 a + (k s d s ) 2 cos 2 a] 
k s £ s [cos 2 a + (k s d s ) 2 sin 2 a] + Vfk s d s 



(21) 
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Fig. 3. T c vs. misorientation angle 2a. The curves cor- 
respond to different thicknesses of the F layers df. The 
parameters are the same as in Fig. 2. 

For the P and AP alignments, under additional as- 
sumption of strong ferromagnetism (h ^> nT cs ), we ob- 
tain: 

^-*G + t&&)-*G). w 

where ip is the digamma function, Vh is determined by 
Eqs. (13) with k\ = (1 + i)kh, and in the region of pa- 
rameters, where T c ^ [the corresponding conditions 
can be extracted from the results for the critical thick- 
ness — see Eqs. (25), (26) below], we may write 

W = RcV h + (ImV h fd s /i s . (24) 

Due to symmetry, the result for the P case (22) repro- 
duces that for the SF bilayer with S layer of thick- 
ness d s [13]. In the AP case, if the second terms 
in the r.h.s. of Eq. (24) can be neglected (e.g., at 
khdf 3> 1 in the region of parameters where T c ^ 0), 
then W = ReVh and we reproduce the result of Ref. 
[8]. However, the second term becomes essential in the 
Cooper limit, defined by conditions d s -C \J D s /2wd, 
dj <C vtun{yj D f / 2lo d , k^ 1 ), 7& = 0, with lod the Debye 
energy of the S material. In this case Re 14 =0 and 
Eqs. (23), (24) reproduce the result of Tagirov [5]. 

The critical thickness d sc of the S layer, below 
which the superconductivity vanishes, immediately fol- 
lows from Eqs. (22), (23) for the P and AP cases: 

df c /6 = 2e c \V h \ , dff/Z. = 2e c W (25) 

at 

d TC /6 < 1. (26) 



Here C « 0.577 is Euler's constant. Condition (26) is 
necessary for applicability of Eqs. (25). If this condi- 
tion is not satisfied, then Eqs. (25) only tell us that 
at dg/^g <C 1 the superconductivity is certainly ab- 
sent, i.e., T c = 0. According to the monotonic growth 
of T c (a), the function d sc (a) decreases monotonically, 
hence df c > d^f. At j b = 0, k h d f > 1, Eqs. (25) re- 
produce the results of Ref. [11] for the P and AP cases. 

The T c (a) dependence can be most easily studied in 
the Cooper limit. In this case a simple analysis (see, 
e.g., Appendix Al in Ref. [13]) can be done already on 
the level of the Usadel equations, and the system is de- 
scribed as a uniform layer with the effective exchange 
energy 2 ^ 

h e g = (T f /T s )h cos a, (27) 

where T s tf) — 2d s (f ) RbA/ 'Ps(f)D s (f)- The accuracy of 
this result is limited to the first order over h, which 
becomes insufficient in the vicinity of a — it/2. At 
a = it/2, the first-order effect of h vanishes, while a 
more accurate analysis (Ref. [5] and Eqs. (23), (24)) 
reveals the second-order effect of h on T c . 

Let us now consider the same limit as in Ref. [11]: 

fe h d/»l, /i>7rT cs , 76 = 0, (28) 

-ykhtfdJZ, « l. (29) 

The condition to have superconductivity at least at 
some orientations has the form df^ < d s <C £ a > an d 
in the case under discussion, Eqs. (25), (26) yield: 

2e c 1 k h £ f < d s /£ s « 1, (30) 

hence condition (29) becomes redundant. 

Starting from Eqs. (20), (15), (21), we finally obtain 
the following equation for T c : 

T c ^ v \2 2nT c ) ^\2 2-kTJ * \2 ) ' 

(31) 

where 

1 • 2 

„ 1 sin a „ 

Q = - + ^^ R=l-Q, 
1 2vsin a — 4 cos 2 a 

£l\.2 = ^-nT cs ^1 + cos 2 a ± V sin 4 a — 4 cos 2 aj , 

do = lk h ^ s /2. (32) 

In the P and AP cases, where the triplet component is 
absent, Eqs. (31), (32) reproduce the results of Rcfs. 
[6,11]. At the same time, at a noncollinear alignment 
the results are clearly different . 

2 ) Since u> n was neglected in comparison with h in the Usadel 
equation, the result of the Cooper limit is valid only at t s ~3> Tf. 
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Fig. 4. Critical thickness of the S layer d sc vs. misorien- 
tation angle 2a. Dashed line is the result of Ref. [11], 
obtained without account of the triplet component. 



The critical thickness is found from Eqs. (31), (32): 
d sc (a)/d = 4V2e c cos ax (33) 

sin 2 a 

(1 + cos 2 a + \J sin 4 a - 4 cos 2 a \ 2 V sin4 Q - 4cos2 Q 
1 + cos 2 a — \J sin 4 a — 4 cos 2 a J 

Although the square root in this expression can become 
imaginary, the whole expression remains real (z l is real 
if \z\ = 1). Figure 4 illustrates the result (33). 

Now we turn to analyze the conditions of applicabil- 
ity for the results reported in Ref. [3]. A noncollincar 
FSF trilayer is a unit cell of the multilayered structure 
studied in that work. The main result of Ref. [3], the 
Josephson current due to the long-range triplet compo- 
nent, requires that the S layer is thin d s <C £ s , while the 
F layers are thick for the singlet component and moder- 
ate for the triplet one: k^ 1 <C £/ < df [3]. In this case 
the condition that superconductivity is not completely 
suppressed at least in the vicinity of the AP alignment 
[Eqs. (25), (26)] takes the form 

4e ^ (TTWTI (34) 

At 7;, = (as it was assumed in Ref. [3]), this yields 
2e c -fkh^f < dg/^s <C 1, which is a rather strong condi- 
tion for 7, since kh£,f ~> 1. Finite interface transparency 
relaxes this condition: already at j b ~ 1, Eq. (34) yields 
2e c 7/ 7fc < d s /i s < 1. 

The condition that superconductivity exists at all 
orientations has the form similar to Eq. (34) but with 
the corresponding expression for df c instead of dff 
in the l.h.s. This only leads to a minor difference, 
since the two critical thicknesses are of the same order: 
df c = \/2df c p at 7 b = 0, while df c = dff at j b > 1. 

In conclusion, we have studied T c of a FSF trilayer 
as a function of its parameters, in particular, the an- 



gle between magnetizations of the F layers. The angu- 
lar dependence becomes pronounced when the S layer 
is thin, and can lead to switching between supercon- 
ducting and non-superconducting states as the angle is 
varied. Our results directly apply to multilayered SF 
structures, where a FSF trilayer is a unit cell. We have 
formulated the conditions which are necessary for exis- 
tence of recently investigated odd triplet superconduc- 
tivity in SF multilayers [3]. 
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